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Abstract 

This article considers some classes of models dealing with the dynamics of discrete 
curves subjected to stochastic deformations. It turns out that the problems of 
interest can be set in terms of interacting exclusion processes, the ultimate goal 
being to derive hydrodynamic limits after proper scalings. A seemingly new method 
is proposed, which relies on the analysis of specific partial differential operators, 
involving variational calculus and functional integration: indeed, the variables are 
the values of some functions at given points, the number of which tends to become 
infinite, which requires the construction of generalized measures. Starting from 
a detailed analysis of the ASEP system on the torus Z/NZ, we claim that the 
arguments a priori work in higher dimensions (abc, multi-type exclusion processes, 
etc), leading to sytems of coupled partial differential equations of Burgers' type. 



Keywords Cauchy problem, exclusion process, functional integration, hy- 
drodynamic limit, martingale, weak solution. 



1 Preliminaries 

Interplay between discrete and continuous description is a recurrent question 
in statistical physics, which in some cases can be answered quite rigorously 
via probabilistic methods. In the context of reaction-diffusion systems, this 
is tantamount to studying fluid or hydrodynamics limits. Number of ap- 
proaches have been proposed, in particular in the framework of exclusion 
processes, see e.g. [T5J,[4| [18], |J3] and references therein. As far as fluid or 
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2 Model definition 



hydro dynamic limits are at stake, most of these methods have in common to 
be limited to systems for which the stationary states are given in closed prod- 
uct forms, or at least for which the invariant measure for finite N (the size of 
the system) is explicitly known. For instance, ASEP with open boundary can 
be described in terms of matrix product form (a sort of a non-commutative 
product form) and the continuous limits can be understood by means of 
brownian bridges (see We propose to address this question from the fol- 
lowing different point of view: starting from discrete sample paths subjected 
to stochastic deformations, the ultimate goal is to understand the nature of 
the limit curves when N increases to infinity. How do these curves evolve 
with time, and which limiting process do they represent as t goes to infinity 
(equilibrium curves)? Following |10| ITT| 1 12 j . we plan to give some partial 
answers to these problems in a series of forthcoming papers. 

The method proposed in the present study is applied in detail to the ASEP 
model. The mathematical kernel relies on the analysis of specific partial 
differential equations involving variational calculus. A usual sequence of 
empirical measures is shown to converge in probability to a deterministic 
measure, which is the unique weak solution of a Cauchy problem. Here 
variables are in fact the values of some function at given points, and their 
number becomes infinite. 

In our opinion, the approach presents some new features, and very likely ex- 
tend to higher dimensions, namely multi-type exclusion processes. A future 
concern will be to establish a complete hierarchy of systems of hydrodynamic 
equations, the study of which should allow us to describe non-Gibbs states. 
All these questions form also the matter of ongoing works. 

2 Model definition 

2.1 A general stochastic clock model 

Consider an oriented sample path of a planar random walk in 1Z 2 , con- 
sisting of N steps (or links) of equal size. Each step can have n discrete 
possible orientations, drawn from the set of angles with some given origin 
{9k = ^-,k = 0, ...,n — 1}. The stochastic dynamics in force consists 
in displacing one single point at a time without breaking the path, while 
keeping all links within the set of admissible orientations. In this operation, 
two links are simultaneously displaced, what constrains quite strongly the 
possible dynamical rules 



2.2 Examples 
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2.1.1 Constructing a related continuous-time Markov chain 

• Jumps are produced by independent exponential events. 

• Periodic boundary conditions will be assumed, this point being not a 
crucial restriction. 

• Dynamical rules are given by a set of reactions between consecutive 
links, an equivalent formulation being possible in terms of random 
grammar. 

With each link is associated a type, i.e. a letter of an alphabet. Hence, for 
any n, we can define which later on will be sometimes referred to as a local 
exchange process. 

For i G [1, N] and k G [1, n\, let Xf represent a link of type k at site i. Then 
we can define the following set of reactions. 

X*X\ +1 +± X\X* +1 , k = l,...,n, l^k + -, 
k Sk+i 

The red equations does exist only for even n, because of the existence of 
folds [two consecutive links with opposite directions], which yield a richer 
dynamics. 

Xf can also be viewed as a binary random variable describing the occupation 
of site i by a letter of type k. Hence, the state space of the system is 
represented by the array 

v ^{X*,i = l,...,N;k = l,...,n}. 



2.2 Examples 

(1) The simple exclusion process 

The first elementary and most sudied example is the simple exclusion pro- 
cess: this model, after mapping particles onto links, corresponds to a one- 
dimensional fluctuating interface. Here we take a binary alphabet and, let- 
ting X 1 = t and X 2 = f , the set of reactions simply rewrites 

A- 

rf ^ fr, 

A+ 

where \^ are the transition rates for the jump of a particle to the right or 
to the left. 



4 



3 Hydrodynamics for the basic asymmetric exclusion process (ASEP) 



(2) The triangular lattice and the ABC model 

Here the evolution of the random walk is restricted to the triangular lattice. 
Each link (or step) of the walk is either 1, e 2i7r//3 or e 4i?r//3 , and quite naturally 
will be said to be of type A, B or C, respectively. This corresponds to the 
so-called ABC model, since there is a coding by means of a 3-letter alphabet. 
The set of transitions (or reactions) is given by 

AB h BA, BC h CB, CA ±=> AC, (2.2) 

p+ q+ r + 

where there is a priori no symmetry, but we will impose periodic boundary 
conditions on the sample paths. This model was first introduced in [8] in 
the context of particles with exclusion, and for some cases corresponding to 
the reversibility of the process, a Gibbs form for the invariant measure was 
given in [9] 

3 Hydrodynamics for the basic asymmetric exclusion process 
(ASEP) 

As mentioned above, we aim at obtaining hydrodynamic equations for a class 
of exclusion models. The method, although relying on classical powerful 
tools (martingales, relative compactness of measures, functional analysis), 
has some new features which should hopefully prove fruitful in other contexts. 
The essence of the approach is in fact contained in the analysis of the popular 
ASEP model, presented below. We note the difficulty to find in the existing 
literature a complete study encompassing various special cases (symmetry, 
total or weak asymmetry, etc). 

Consider iV sites labelled from 1 to N , forming a discrete closed curve in the 
plane, so that the numbering of sites is implicitly taken modulo N , i.e. on 
the discrete torus G (JV) = Z/AZ. In higher dimension, say on the lattice Z k , 
the related set of sites would be drawn on the torus (Z/AZ) fc . 

We gather below some notational material valid throughout this paper. 

• 1Z (resp. 1Z + ) stands for the real (resp. positive real) line. C k [0, 1] is 
the collection of all real-valued, fc-continuously differentiable functions 
defined on the interval [0, 1], and Ai is the space of all finite positive 
measures on the torus G = [0, 1). 

• For S an arbitrary metric space, V(S) is the set of probability mea- 
sures on S, and T>s[0,T] is the space of right continuous functions 
z : [0, oo] — > S with left limits and t — > z t . 
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• C^{K) is the space of infinitely differentiable functions with compact 
support included in K £ S, and we shall write C[T] to denote the 
subset of functions <fi(x,t) € Cq°([0, 1] x [0, T_]) vanishing at t = T. 

• For i = 1, . . . , N, let A^it) and B!- N '(t) be binary random variables 
representing respectively a particule or a hole at site i, so that, owing 
to the exclusion constraint, A[ (t) + B^ N \t) = 1, for all 1 < i < iV. 
Thus (AW(i) = (Af . . . , Aj^(t)),t > 0} is a Markov process. 

• f2 (JV) will denote the generator of the Markov process A (JV) (i), and 
J-f: = (t(A (jv) (s), s < i) is the associated natural filtration. 

• Our purpose is to analyze the sequence of empirical random measures 

Mf 5 = ^ E (3.1) 

when A — > oo, after a convenient scaling of the parameters of the 
generator f2 (JV) . The probability distribution associated with the path 
of the Markov process [if ,t € [0, T], for some fixed T, will be simply 
denoted by Q<M. 

As usual, one can embed G (JV) in G, so that a point i € G (JV) corresponds 
to the point i/N in G. Hence, in view of (|3.ip . it is quite natural to let 
the sequence Q (JV) be defined on a unique space Dm[0,T], which becomes a 
polish space (i.e. complete and separable) via the usual Skorokod topology, 
as soon as Ai is itself Polish (see e.g. [7j, chapter 3). Without further 
comment, Ai is assumed to be endowed with the vague product topology, as 
a consequence of the famous Banach-Alaoglo and Tychonoff theorems (see 

e.g. nans]). 

Choose two arbitrary functions (j) a ,(j>b G C[T] and define the following real- 
valued positive measure 



Zi N> [4> a ,4>b] = exp 



(3.2) 



viewed as a functional of (j) a , <f> h . Since Af\t) + B\ N) {t) = 1, for 1 < i < N, 
the transform Z\ ' is essentially a functional of the sole function (j) a — (f)^, 
up to a constant uniformly bounded in N. Nevertheless, it will appear later 
that we need 2 independent functions. For the sake of brevity, the explicit 
dependence on N,t or <fi, of quantities like A[ N \t), B^ N \t), Z t (JV) [cj) a , (fib], will 
frequently be omitted, wherever the meaning remains clear from the context: 
for instance, we simply write Ai,Bi or Z[ N \ Also Z (JV) stands for the process 

{zr,t>oy 

A standard powerful method to prove the convergence (in a sense to be 
specified later) of the sequence of probability measures introduced in (|3.ip 
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3 Hydrodynamics for the basic asymmetric exclusion process (ASEP) 



consists first in showing its relative compactness, and then in verifying the 
coincidence of all possible limit points (see e.g. [14] . Moreover here, by the 
choice of the functions <t> a ,<l>bi it suffices to prove these two properties for 
the sequence of projected measures defined on T>-ji[0,T] and corresponding 
to the processes {Z[ N [4> a , fib], t > 0}. 

Let us now introduce quantities which, as far as scaling is concerned, are 
crucial in order to obtain meaningful hydrodynamic equations. 

" def A a6 (JV) + A 6 °(iV) 
MN) ~ ~ 2 > (3>3) 

(i(N) = X ab (N)- X ba (N), 

where the dependence of the rates on N is explicitly mentioned. 

Theorem 3.1. Let the system \3. 3\) have a given asymptotic expansion of the 
form, for large N, 

'\(N) = XN 2 + o(N 2 ), 

(3.4) 

H(N) = fiN + o(N), 

where X and p are fixed constants. [As for the scaling assumption ( fff.^l j, the 
random measure log Z+ is a functional of the underlying Markov process, 
in which the time has been speeded up by a factor N 2 and the space shrunk 
by N^ 1 ]. Assume also the sequence of initial empirical measures log^Q , 
taken at time t = 0, converges in probability to some deterministic measure 
with a given density p(x,0), so that, in probability, 



limlogZ^/ [p(x,O)0 a (x,O) + (l-p(x,O))0 b (x,O)]dx, (3.5) 



for any pair of functions (j) a , (fib £ C[T]. 

Then, for every t > 0, the sequence of random measures converges in 
probability, as N — )• oo, to a deterministic measure having a density p(x,t) 
with respect to the Lebesgue measure, which is the unique weak solution of 
the Cauchy problem 
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Ot 



dx 2 



Hp(x,t)(l - p{x,t)) 



d<fi(x, t) 
dx 



dxdt 



+ / p(x,0)(j)(x,0)dx = 0, 



(3.6) 



where \3. b)) holds for any function (f) G C[T]. 



3.1 Existence of limit points: sequential compactness 
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If, moreover, one assumes the existence of d g^i°\ f or p( x >fy given, then 
\3. 6\) reduces to a classical Burgers' equation 

dp(x,t) x d 2 p{x,t) dp(x,t) 
—&T~ = dx? + ^ ~ 2P{X > t)] ^x— 

Proof. The proof is contained in the next three subsections. 



3.1 Existence of limit points: sequential compactness 



As usual in problems dealing with convergence of sequences of probability 
measures, our very starting point will be to establish the weak relative com- 
pactness of the set {log Z[ N \ N > 1}. Some of the probabilistic arguments 
employed in this paragraph are classical and can be found in the literature, 
e.g. |18[ 114] . although for slightly different or simpler models. 

Letting 4> a ,4>b be two arbitrary functions in C[T], we refer to equation (|3.2p . 

Using the exponential form of Z t (JV) and Lemma [Al-5.1] in |14j (see also chap- 
ter 3 in [7] for related calculus), one can easily check that the two following 
random processes 

jjw def z w_ z m_ I (n^[z ( s N) ] + ei N) z^)ds, (3.7) 

J o 



(iV) def ( TT (N)\2 



(V? 



flWpW) 2 ] - 2ZWnW[ZW])da (3.8) 



are bounded {J-} }-mar ting ales, where 



o(JV) def 

* ~ N 



dt \N 



(3.9) 



Setting now 



ib = 

Yxy — 



we have 



CXT>( — Alb. xy y— . 



where 



xy = ab or ba, 



(3.10) 



L ( t N) = E ^ b ) (^)^B t+1 + XZ\i,t)B t A l+1 . (3.11) 
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3 Hydrodynamics for the basic asymmetric exclusion process (ASEP) 



On the other hand, a straightforward calculation in equation (|3.1ip allows 
to rewrite (|3.8p in the form 

V W = ([/ W )2 _ h Z WfR^)ds, (3.12) 
Jo 

where the process is stricly positive and given by 

The integral term in fj3. 12[) is nothing else but the increasing process associ- 
ated with Doob's decomposition of the submartingale (C/ t (JV) ) 2 . 

The folllowing estimates are crucial. 

Lemma 3.2. 

L[ N) = 0(1), (3.13) 

R\ N) = O(^). (3.14) 

Proof. We will derive fj3. 13|) by estimating the right-hand side member of 
equation (|3.1ip . From now on, for the sake of shortness, the first and second 
partial derivatives of tp(z, t) with respect to z will be denoted respectively 
by ip'(z,t) and ip"(z,t). 

Clearly, Aip xy ^,t) = jj^xyijf^) + ^(jv 5 )' Then, taking a second order 
expansion of the exponential function and using equations fj3.3[) and (|3.4p . 
we can rewrite (|3. as 



T (n) _ K N ) sr 







+ E (Ai-A i+1 )A^ ab ^,t)+0^). (3.15) 



The first sum on the right in (I3.15P is uniformly bounded by a constant 
depending on tp. Indeed \Ai\ < 1, and ip' is of bounded variation since 
ip 6 C[T]]. As for the second sum coming in (|3.15p . we have 

Then the discrete Laplacian 



AVa6(^-,*) -A^ a6 (-^,i) 



3.1 Existence of limit points: sequential compactness 
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admits of the simple expansion 

A^(i±i,t) - A^(±,t) = ^{p) +0(±). (3.16) 

By ([321), K N ) = xn2 + o(N 2 ), so that fl3~T6l) implies 



= 0(1), (3.17) 

which concludes the proof of (|3.13p . The computation of R^' leading to 
f!3.14p can be obtained via similar arguments, remarking that 

Rf^aAb] = 412^,24] - 2L\ N) [cp a ,cP b ]. 



To show the relative compactness of the family Z {N \ which from the sepa- 
rability and the completeness of the underlying spaces is here equivalent to 
tightness, we proceed as in |14] by means of the following useful criterion. 



Proposition 3.3 (Aldous's tightness criterion, see [1J). A sequence {X (N) } of 
random elements ofT>n[0,T] is tight (i.e. the distributions of the {X (JV) } are 
tight) if the two following conditions hold: 



(i) 



lim lim sup P[\ \X (N) \ \ > a] = 0, (3.18) 

(2— >-00 



where \\X™\\ = sup |X t (JV) 1 



t<T 



(ii) For each e,7], there positive numbers 5q and Nq, such that, if 5 < 5o 
and N > Nq, and if r is an arbitrary stopping time with t + 6 < T, 
then 

P[\X^-XW\>e] < V . (3.19) 

Note that condition \'J.18\) is always necessary for tightness. ■ 

We shall now apply Lemma 13.21 to equations (|3.7p and (|3.12p , the role of 
Xj. in Proposition 13.31 being played by Z[ N) . 

The random variables Z 4 (JV) and 9^ are clearly uniformly bounded, so that 
condition (|3.18p is immediate. To check condition (|3.19p . rewrite (|3.7p as 



~ A t - u t+8 



U { t N) + J* +& (Lf) + 9^)Z^ds. (3.20) 
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3 Hydrodynamics for the basic asymmetric exclusion process (ASEP) 



For t £ [0, T], the integral term in (|3.20p is bounded in modulus by Kb 
{K being a constant uniformly bounded in N and ifi), and it satisfies (|3.19p . 
whenever t is replaced by an arbitrary stopping time. We are left with the 
analysis of U[ . But, from (|3.12p . (|3.1 4p and Doob's inequality for sub- 
martingales, we have 



E[«-^) 2 ] 



t+s 



Zi N) ) 2 R[ N) ds 



sup \Ut I > e 

t<T 



< 4e 



Jo 



< 



< 



Cb_ 

Jf' 

ACT 



(3.21) 



where C is a positive constant depending only on ip. Thus — > in 
probability, as N — > oo. This last property, together with (|3.7p . (|3.20p and 
assumption (|3.5p . yield (|3.19p and the announced (weak) relative compact- 
ness of the sequence Z[ N) . Hence, the sequence of probability measures Q (N \ 
defined on Pvi[0,T] and corresponding to the process fif N \ is also relatively 
compact: this is a consequence of classical projection theorems (see for in- 
stance Theorem 16.27 in |13]). We are now in a position to establish a further 
important property. 



Let Q the limit point of some arbitrary subsequence Q^ Nk \ as — > oo, 
and Zt = limjVj._j.oo z[ Nk \ Then the support of Q is a set of sample paths 
absolutely continuous with respect to the Lebesgue measure. Indeed, the 
application fa — > sup t< <p log Zt is continuous and we have the immediate 
bound 



sup log Z t < 

t<T 



t>a(x,t)\ + \4> b (x,t)\]dx, 



which holds for all ipa^ipb £ C 2 [0, 1]. Hence, by weak convergence, any limit 
point Zt has the form 



Zt[4>a,<t>b] = exp / [p(x,t)4> a (x,t) + (1 - p(x,t)(p b (x,t)]dx 
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(3.22) 



where p(x,t) denotes the limit density [which a priori is a random quantity] 
of the sequence of empirical measures pi introduced in (|3.ip . 



3.2 A functional integral operator to characterize limit points 

This is somehow the Gordian knot of the problem. Relying on the above 
weak compactness property, our next result shows that any arbitrary limit 
point Q is concentrated on a set of trajectories which are weak solutions of 
an integral equation. 



3.2 A functional integral operator to characterize limit points 
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The main idea is to consider for a while the 2N quantities cj) a {jj, t) , (j)b{jj,t) , 
1 < i < N, as ordinary free variables, which for the sake of shortness will be 
denoted respectively by x[ and y\ N ^ ■ With this approach, the problem of 
the hydrodynamic limit will appear to be mostly of an analytical nature. 



Let 



<>(M) = A ab (iV) 



<>(M) = A ba (iV) 



cxp 



exp 



x l+i ~ x T J + v\ 



N 

(N) _ „,(JV) „(N) _ (JV) 
X j+1 



Vi+i ~Vi + x 



N 



1 



Then, using (|577]l . (ETUI) . (l3~TU1) . (I37TT1) and the definition of Z 4 (iV) , we obtain 
immediately the following functional partial differential equation (FPDE) 

C[ N) [Z ( t N) } + ei N) Z { t N \ (3.23) 



d(zr-un 



dt 

where is the operator 



d 2 h 



More precisely, introducing the family of cylinder sets 

V W = [-|*U*I] P . P = l,2..., (3.24) 



with 



|*|= sup (\(/) a (z,t)\,\(f) b (z,t)\), 

(z,t)e[o,i]x[o,T] 

,(JV) 



(3.25) 



we see immediately that C\ acts on a subspace of CQ°(V (2Ar) ), since Z\ is 
analytic with respect to the coordinates {0 O (., i), i)}, for each finite ./V 
(things will be made more precise in section I3.2.3[) . Needless to say that 
is not of parabolic type, as the quadratic form associated with the second 
order derivative terms is clearly non definite (see e.g. [6]). In addition, equa- 
tion (|3.23p is a well defined stochastic FPDE, as all underlying probability 
spaces emanate from families of interacting Poisson processes. 

Now, it might be worth accounting for the use of the word functional above, 
and for the reason of isolating the third term on the right in (|3.23p . Indeed, 
by (|3.9p . Oj. is a functional involving also a partial derivative of Z[ N) with 
respect to t, as we can write 



(JV) 7 (N) 



i 



- y 



\dZ, 



(JV) 



N ^ ldx (N) dt 



dx^> + dz { t "> d y . 



(JV) 



dyl N) dt 



dZ\ 



(JV) 



dt 



(3.26) 
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3 Hydrodynamics for the basic asymmetric exclusion process (ASEP) 



The last equality in (|3.26p might look somewhat formal, but will come out 
more clearly in Section [3.2,31 Note also, by the same argument which led to 
(|3.22p . that we have 

lim 6 { t Nk) = [\p( x ,t)^(x,t) + (l-p(x,t)^(x,t)]dx. (3.27) 
iVfc— >oo Jq ot at 



Our essential agendum is to prove that any limit point of the sequence of ran- 
dom measures pt w = k li m p{ Nk ' satisfies an integral equation correspond- 

ing to a weak solution (or distributional in Schwartz sense) of a Cauchy type 
operator. To overcome the chief difficulty, namely the behaviour of the limit 
sum in (|3.23p . we propose a seemingly new approach, which is tantamount 
to analyzing the family of second order linear partial differential operators 
£j JV) along the sequence JVjt — > oo. 

As briefly emphasized in the remark at the end of this section, a brute force 
analysis of (|3.23p would lead to a dead end. Indeed an important preliminary 
step consists in extracting the juice of the estimates obtained in Lemma [3.21 
to rewrite the operator C$ in terms of only N principal variables, up to 

quantities of order O^j^j. This is summarized in the next lemma. 
Lemma 3.4. The following FPDE holds. 

a\A, t u t ) def , m , fl (») 7 (jv) ,m( L 



(U Af> [Z ^]+er>Zr + 0[^), (3.28) 

where A[ N) is viewed as an operator with domain Cq°(V N ) such that 



1 / dg dg \ d 2 g 



-I f_ — TV- 

2 \dx[ N) dx[1\J dxf^dx^-y 

Ofj 



(3.29) 



the term O(^) being in modulus uniformly bounded by jj, where C denotes 
a constant depending only on the quantity 

sup {ip(x,t),ip'(x,t),ip"(x,t)} . 

(x,t)€[0,l]x[0,T] 

Proof. The result follows by elementary algebraic manipulations from equa- 
tions §3]), (133), fl3~TUj) . pT31) . pTfjl . and details will be omitted. ■ 



3.2 A functional integral operator to characterize limit points 
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Taking Lemma [3.41 as a starting point, we sketch out below in Sections 13.2.11 
and 13.2.21 the main lines of our analytical approach, which indeed can be 
briefly summarized by means of some keywords. 

• Coupling, taken here in Skohorod's context. 

• Regularization and functional integration. Regularization refers to the 
fundamental method used in the theory of distributions to approxi- 
mate either functions or distributions. We shall apply it to the FPDE 
(|3.28j) . considering t and the N values of the function O (., t) (taken at 
points of the torus) as N + 1 ordinary variables. Then, passing to the 
limit as N — > oo, we introduce convenient functional integrals together 
with variational derivatives. This might likely extend to much wider 
systems, although this assertion could certainly be debated. 



3.2.1 Interim reduction to an almost sure convergence setting 

This can be achieved by means of the extended Skohorod coupling (or trans- 
fer) theorem (see Corollary 6.12 in [13J), which in brief says that, if a sequence 
of real random variables is such that lim^oo fk(£,k) = /(£) converges 
in distribution, then there exist a probability space V and a new random 
sequence such that = and lim^oo fk(Ck) = /(£) 5 almost surely 
in V, with £ = £. Here this theorem will be applied to the family z[ Nk \ 
which thus gives rise to a new sequence, named Y^ Nk ' in the sequel. Clearly 
this step is in no way obligatory, but rather a matter of taste. Indeed, one 
could still keep on working in a weak convergence context, with Alexandrov's 
portmanteau theorem (see e.g. [7]) whenever needed. ■ 

3.2.2 Considering p.290 as a partial differential operator with 
constant coefficients 

For each finite N, we consider the quantities 




as constant parameters, while the x\ 's will be taken as free variables from 
a variational calculus point of view. This is clearly feasible, remembering 
that, by definition, ip a b = 4> a — <ftb, for all 4> a ,4>b £ C[T]. Hereafter, t will be 
viewed as an exogeneous mute variable, not participating concretely in the 
proposed variational approach. 

Then, according to the notation introduced in section 13.2.11 we can rewrite 
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(l3T28j) in the form 

d[Yt - Ut ' d = f A ( t N) [Y t (N) ) + el N) Y t (N) +°{jf)> ( 3 - 3 °) 

where 6[ N) is still given by (|3.9p , keeping in mind that all random variables in 
(|3.30p are defined with respect to this new (though unspecified) probability 
space introduced in section 13.2,11 above. In particular, from the tightness 
proved in section 13.11 letting N — > oo along some subsequence iV& , we have 

lim r t (iVfe) [0 a! 6 ] a 4-F t [0 o ,0 6 ]. 

fe— >oo 



3.2.3 Analysis of the FDPE flX3T)ft 

The idea now is to propose a regularization procedure, which consists in 
carrying out the convolution of (|3.30p with a suitably chosen test function, 
noting that Y^ nowhere vanishes and is uniformly bounded. 

As usual, the convolution / * g of two integrable functions /,j G Co°(V (Ar) ), 
with V (JV) given in (|3.24p . will be defined by 



{f*g){u)= f(u-v)g(v)dv. 
Let to be the function of the real variable z defined by 

exp(~) if z < 0, 



if z > 0, 



def dw(z) , II , s def d 2 Uj(z) 



where it will be convenient to write U)'(z) = £ z and uj"{z) — -7-3 

Setting f<"» = {x[ N \x ( 2 N) ,...,x^ } ), with x\ N) = a (£,t),l < % < N, we 
introduce the following family of positive test functions Xe^ £ Co°(V (JV) ), 

xfH^ = ^X>r) 2 -A e>0. (3.31) 
^ i=l ' 



Note: It is worth keeping in mind that, as often as possible, the time vari- 
able t will be omitted in most of the mathematical quantities, e.g. x (JV) (i). 
Indeed, as mentioned before, t plays in some sense the role of a parameter. 

From (|3.30p . it follows immediately that 

' d(Y {N) - n (N) ) \ 

(3.32) 
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Now, by (|3.29p and according to the statement made in Section 13.2.21 the 
first term in the right-hand side member of f|3.32[) can be integrated by parts, 
so that, for any finite N and e sufficiently small, 

(A N) {Yl N) ]*xi N} ) (x (JV) ) = (A N) lxr]*Yt iN) ) (* (A °), (3-33) 

where A^ is by definition the adjoint operator of .A (JV) in the Lagrange 
sense. Here the domain of Af consists of all functions /i (iV) of the form 



h w 



exp 



dai N) (x)V[<f> a (x, t)\ + da { , N) (x)V[<f> b (x, t)\ 



where 



V : C[T] — > 1Z + stands for an arbitrary analytic function; 

and crJ^' are arbitrary discrete probability measures on G (iV) . 



Clearly Xe belongs to the domain of A t . Under the assumptions made in 
Section [3.2,2[ a direct integration by parts in (|3.29p yields the formula 



1 



2 \dx 



dh dh 
+ 



(N) 



ux i+l 



+ N 



d 2 h 



dx\ N) dx^ 



dh 

A >d^ 



(3.34) 



where A ( t N ' > has been defined in (|3.33p . 



Let, for each 



Xe 



b a E C[T], 



lim xr(x {N) )=u 



4> 2 {x,t)dx — e' 



(3.35) 



where the integral in f|3. 35[) is readily obtained as the limit of the Riemann 
sum in (|3.3ip . 

Lemma 3.5. For each 4>(x,t) € C[T], the following limit holds uniformly. 



lim Al N) [xi N) }(^) 
N^foo 



[^'abi x i t) K (4>, X, t) + \tpab(x, t)H((j), X, t)] dx, 

(3.36) 



with 



H((f>, z, t) = 2<f>(z, t)J ( J <f> 2 (u, t)du - e 2 ^j . 



K(<j), z, t) = H((j), z, t) + 4</> 2 (z, t)uj' : 



(j) 2 (u,t)du — e 2 ) . 



(3.37) 
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3 Hydrodynamics for the basic asymmetric exclusion process (ASEP) 



Proof. For fixed N, the function -A^fXe^] has partial derivatives of any 
order with respect to the coordinates x\ , i = 1,... ,N. In particular, we 
get from (ROB , for each i G G (JV) , 

?)V {N) T (N) f 1 N n \ 

- 2 V'(^i>n 2 -4 

! V 1=1 7 



fl2 v (iV) T W (AO / -, JV x 



v i=l 



So we are again left with Riemann sums, which, for any continuous function 
k(x,t), yield at once 

d\i N) 



[X> 



j'(^J (/) 2 (x,t)d. 



l 

2 



2u/ / </r(a;,i)dx-ir / k(x,t)c/)(x,t)dx 



o 



and 



^"(J 4> 2 (x,t)dx - e 2 ^ J k(x,t)(f> 2 (x,t)dx. 

Keeping in mind that, for all N, the vector x (JV) (from its very definition) 
must be a discretization of some function element in C[T], the last important 
step to derive (|3.36p requires to give a precise meaning to the limit 



lim / x ( £ N \^ Nr )dx m , 
N-XxJ v (n) 

allowing to carry out functional integration and variational differentiation. 
In this respect, let us emphasize (if necessary at all !) that the usual con- 
structions of measures and integrals do not apply in general when the domain 
of integration is an infinite-dimensional space of functions or mappings, all 
the more because a complete axiomatic for functional integration does not 
really exist; indeed each case requires the construction of ad hoc generalized 
measures, see promeasures in [21 [3] or quasi-measures in |16J , These ques- 
tions prove to be of importance in various problems related to theoretical 
physics. In our case study, integration over paths belonging to G C[T] needs 
to be properly constructed. 

The point is to define, as N — > oo, a volume element denoted by 5 (</>). 
We shall not do it here, but this could be achieved by mimicking classical 
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fundamental approaches, see e.g. [21 El [16]. The main tool is the important 
F. Riesz' representation theorem, which to every positive linear functional 
let correspond a unique positive measure. ■ 

Now from equation fj3. 35[) it is permissible to introduce the normalized test 
functional 

- f*\ Xe(0) 
XeW) = —fi-i 

where D is chosen to ensure 

cg°([-l*l,l*|]) 
and is given by (|3.25p . so that 

D= f co ( f tf(x,t)dx-£ 2 \ 5{(f>). 

7c§°([-|*|,|*|]) V^o / 

From Skohorod's coupling theorem, Yt does satisfy an equation of the form 
(|3.22p . Hence, we can write the following functional derivatives (which are 
plainly of a Radon-Nykodym nature) 



,2 



p\,t)Y t . 



Now everything is in order to complete the puzzle, according to the following 
steps. 

1. First, using (|3.33p . rewrite (|3.32p as 



+((eryn*xi N w N) )+o(± 



(3.39) 



2. Let N — > oo in (|3.39p and then replace Xei^) by Xe(4>), remembering 
that by (13T2TT) U { t N) = 0(1/ N) -»• uniformly. 

3. Carry out two functional integration by parts in equation (|3.36p by 
making use of (|3.38p . 

4. Finally, integrate on [0, T], let e — > and switch back to the original 
probability space, where Z[ N) [(f) a , fa], by section [3.2. 1[ converges in 
distribution to Zt[(f) a , (f>b\- this yields exactly the announced Cauchy 
problem (|3.6p . 
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Hence, the family of random measures converges in distribution to a 
deterministic measure m, which in this peculiar case implies also convergence 
in probability. 



3.3 Uniqueness 



The problem of uniqueness of weak solutions of the Cauchy problem (|3.6p 
for nonlinear equation is in fact already solved in the literature. For a wide 
bibliography on the subject, we refer the reader for instance to [6]. The proof 
of Theorem 13.11 is concluded ■ 



4 Conjecture for the n-species model 



We will state a conjecture about hydrodynamic equations for the n-species 
model, briefly introduced in section 12.11 in the so-called equidiffusive case, 
precisely described hereafter. 

Definition 4.1. The n-species system is said to be equidiffusive whenever 
there exists a constant A, such that, for all pairs (k,l), 

N^oo N 2 



Then, letting 



o 



kl 



lim N log 



X kl (N) 
X lk (N) 



we assert the following hydrodynamic system holds. 
dpk 



dt 



A 



d 2 p k + ^_( 
dx 2 dx V 



jk 



+ — I X! a ""PkPl 



l^k 



, k = 1, . . . , n. 



The idea is to apply the functional approach presented in this paper: this is 
the subject matter of an ongoing work. 
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